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The poor convergence of quantum field theory at finite temperature has been one of the main 
obstacles in the practical applications of thermal QCD for decades. Here we briefly review the 
progress of hard-thermal-loop perturbation theory (HTLpt) in reorganizing the perturbative expan- 
sion in order to improve the convergence. The quantum mechanical anharmonic oscillator is used as 
a simple example to show the breakdown of weak-coupling expansion, and variational perturbation 
theory is introduced as an effective resummation scheme for divergent weak-coupling expansions. 
We discuss HTLpt thermodynamic calculations for QED, pure-glue QCD, and QCD with Nf = 3 
up to three-loop order. The results suggest that HTLpt provides a systematic framework that can 
be used to calculate both static and dynamic quantities for temperatures relevant at LHC. 
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I. INTRODUCTION 



Since the beginning of experiments in 1999, the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National 
Laboratory has been doing intensive studies at initial temperatures up to twice the critical temperature for decon- 
finementQ T c ~ 170 MeV. This translates to a strong coupling constant of a s (n — 2-k x 170 McV) = g 2 /(47r) ~ 0.4, or 
cquivalcntly g ~ 2, where /x is the renormalization scale and it relates to the temperature by fj, — 2ttT. The upcoming 
experiments at the Large Hadron Collider (LHC) at European Organization for Nuclear Research are expected to 
yield initial temperatures of 4 — 6 T c , driving the running coupling further down. For the RHIC and LHC experiments 
to have the greatest possible impact on science, it is essential to make as close a connection to the fundamental 
theory of quantum chromodynamics (QCD) as possible. There is an urgent need for theoretical analysis that is based 
rigorously on QCD but which can also make contact with more phenomcnological approaches, particularly in the area 
of thermodynamics and real-time dynamics of QCD at intermediate coupling, g ~ 2. 

We have to be extremely careful when dealing with this intermediately coupled region. Naively, g ~ 2 seems to 
suggest the breakdown of perturbation theory in this region. This is also in line with the observations from the early 
RHIC data that the state of matter created there behaved more like a strongly coupled fluid than a weakly coupled 
plasma [HI]. As a result, the term "quark-gluon plasma" might need to be modified to "quark-gluon liquid", and 
a description in terms of strong coupling formalisms, such as hydrodynamics or AdS/CFT correspondence, might be 
more appropriate. However on the other hand, g ~ 2 is not huge especially when considering that a s = g 2 / (47r) ~ 0.4 
is still a small number. So people have not yet totally lost faith on perturbation theory and as a payback observables 
like jet quenching 0, Q and elliptic flow @ have been able to be described using a perturbative formalism. Therefore it 
seems that a complete understanding of QGP would require knowledge from both strong-coupling and weak-coupling 
formalisms, and in this paper we would focus on the latter approach. 

Thermodynamics describes the bulk properties of matter in or near equilibrium which are theoretically clean and 
well defined. The calculation of thermodynamic functions for finite temperature field theories has a long history. In the 
early 1990s the free energy was calculated to order g 4 for massless scalar </> 4 theory [13, EH, quantum electrodynamics 



Note that the deconfinement transition is actually an analytic crossover [lj and T c represents a temperature around which the thermo- 
dynamic quantities change quickly. 
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FIG. 1. Weak-coupling expansion for the scaled QCD pressure with Nf = 3. Shaded bands show the result of varying the 
renormalization scale /x by a factor of 2 around [i = 2irT. 



(Q ED) (l2l ] and QCD [TiIlLS], respectively. The corresponding calculations to order g 5 were obtained soon afterwards 
jl3T - l20l ]. Recent results have extended the calculation of the QCD free energy by determining the coefficient of the 
g log g contribution [2l[ • For massless scalar theories the perturbative free energy is now known to order g 6 22 1 and 
<? 8 log<7 0. 

Unfortunately, for all the above mentioned theories the resulting weak-coupling approximations, truncated order- 
by-order in the coupling constant, arc poorly convergent and show large dependence on the renormalization scale 
unless the coupling constant is tiny which corresponds to astronomically high temperatures. In Fig. [TJ we show the 
weak-coupling expansion for the QCD pressure with Nf = 3 normalized to that of an ideal gas through order g 5 . The 
various approximations oscillate wildly and show no signs of convergence in the temperature range shown which is 
probed in the ongoing experiments. The bands are obtained by varying the renormalization scale fi by a factor of 2 
around the value fi = 2ttT and we use three-loop running for a s [24j with Ajjg(iV/ = 3) = 344 MeV [251 . In Fig. [2] 
we show the weak-coupling expansion for the QED pressure with Nf = 1 normalized to that of an ideal gas through 
order e 5 , and we see clearly the same poor convergence pattern as the QCD case. Therefore this oscillating behavior 
is not specific to QCD, but a generic feature for hot field theories which actually has been also observed in scalar 
theories [HI EH [lj, [HI, HH, [23[ ■ Due to this subtlety, a straightforward perturbative expansion in powers of a s for 
QCD does not seem to be of any quantitative use even at temperatures many orders of magnitude higher than those 
achievable in heavy-ion collisions. 

The poor convergence of finite-temperature perturbative expansions of thermodynamic functions stems from the fact 
that at high temperature the classical solution is not described by massless particle states. Instead one must include 
plasma effects such as the screening of electric fields and Landau damping via a self-consistent hard-thermal-loop 
(HTL) resummation [26j . The inclusion of plasma effects can be achieved by reorganizing perturbation theory. There 
are several ways of systematically reorganizing the finite-temperature perturbative expansion [27l - l29| . In this paper 
we will focus on the hard-thermal-loop perturbation theory (HTLpt) method (30l - [4ll | . The HTLpt method is inspired 
by variational perturbation theory (VPT) [42TI49I] . HTLpt is a gauge- invariant extension of screened perturbation 
theory (SPT) [50l454| . which is a perturbative reorganization for finite-temperature massless scalar field theory. In 
the SPT approach, one introduces a single variational parameter which has a simple interpretation as a thermal 
mass. In SPT a mass term is added to and subtracted from the scalar Lagrangian, with the added piece kept as 
part of the free Lagrangian and the subtracted piece associated with the interactions. The mass parameter is then 
required to satisfy either a variational or perturbative prescription. This naturally leads to the idea that one could 
apply a similar technique to gauge theories by adding and subtracting a mass in the Lagrangian. However, in gauge 
theories, one cannot simply add and subtract a local mass term since this would violate gauge invariance. Instead, 
one adds and subtracts an HTL effective action which modifies the propagators and vertices selfconsistcntly so that 
the reorganization is manifestly gauge invariant (55| . 
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FIG. 2. Weak-coupling expansion for the scaled QED pressure with Nf = 1. Shaded bands show the result of varying the 
renormalization scale /z by a factor of 2 around /i = 27rT. 

This paper focuses on the resummation of thermodynamics for gauge theories. In Sec. [IT] we use a one-dimensional 
quantum mechanical example, namely the anharmonic oscillator, to demonstrate the breakdown of weak-coupling 
expansion in a simple setup. Then in Sec. IIIII we introduce variational perturbation theory which is an effective 
resummation scheme for divergent weak-coupling expansions. Hard-thermal-loop perturbation theory is introduced 
in Sec. IIV1 where we discuss its formalism and results of thermodynamic calculations to three-loop order for QED, 
pure-glue QCD, and QCD with Nf = 3. We conclude in Sec. |V]with a brief outlook for the real-time and other 
applications of HTLpt. 

II. ANHARMONIC OSCILLATOR 

Before entering the complicated field theory, let us make our life simple by first going back to quantum mechanics. 
In this section we will use an example of one-dimensional anharmonic oscillator to show that weak-coupling expansion 
totally breaks down even for this simple system no matter how small the coupling is, in contrast to daily intuition. 

Now let us consider the ground-state energy E of a simple anharmonic oscillator with potential 

v( X ) = i«v + 9 -/ , (i) 

with u and g some arbitrary positive coefficients. The weak-coupling expansion of the ground-state energy has been 
calculated to all orders by Bender and Wu (5rj| - [59| and the result reads 

oo 

^)=-E^ W (£a) > ( 2 ) 

71=0 

where c^ w are rational coefficients obtainable to all orders 

BW _ J 1 3 21 333 30885 

In the limit n — > oo, c^ w becomes 

hmc« w = (-ir+ 1 7j3>-i)!. (4) 

Due to the factorial growth in Eq. (|4]), the weak-coupling expansion ground-state energy in Eq. ([2J is an asymptotic 
series with zero radius of convergence. The resulting weak-coupling expansion approximations to the ground-state 
energy are plotted in Fig. [3] The vertical axis is the anharmonic oscillator ground-state energy scaled by the cor- 
responding harmonic oscillator ground-state energy E(g)/ E(0), and the horizontal axis is the coupling strength g. 
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FIG. 3. Weak-coupling expansion results of the ground-state energy of an anharmonic oscillator. 



Curves in different colors correspond to the ground-state energy up to various orders in the weak-coupling expansion, 
fnstcad of converging, Fig. [3] shows clearly that the weak-coupling results are oscillating. If keep adding higher order 
terms, the results would be bended up and down more steeply. Finally as n — > oo, the resulting curve would blow to 
infinity right from the origin which demonstrates the meaning of zero radius of convergence. This result is striking 
and highly nontrivial: Intuitively "a small coupling expansion" sounds like "a perturbative expansion" , however the 
above result actually indicates that "small coupling" may not be "perturbative" which is totally against our intuition! 
The anharmonic oscillator therefore provides a simple example for the breakdown of weak-coupling expansion, and 
calls the need of rcsummation. 



III. VARIATIONAL PERTURBATION THEORY 



In order to improve the convergence of the weak-coupling expansion and with the inspiration from the Feynman- 
Kleinert variational approach to path integral fsol ] strong-coupling expansion of variational perturbation theory was 
developed in the 1990s gJSS 

The basic idea of VPT is rather simple: First, the harmonic term of the potential is split into a new harmonic term 
with a trial frequency fi and a reminder: 

uj 2 -> o 2 + (lu 2 - n 2 ) . (5) 

Then the anharmonic potential is rewritten into 

V(x) = ^-x 2 + V int (x), (6) 

with an interaction 

V int (x) = ^(rx 2 +x i ), (7) 

where 

r =-{uj 2 - n 2 ) . (8) 
9 

After this rewriting, a perturbation expansion is carried out at fixed r which generates the new ground-state energy: 

JV 

E N (g,r) = nJ2 c ^)(^y - (9) 

n=0 



2 See also the independent development by Giachetti and Tognetti in Refs. |61|463[I 

3 For the development and other variants of VPT, we refer the reader to Refs. I64l4ll4ll for a far from complete list of early references. 
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FIG. 4. Variational perturbation theory results of the ground-state energy of an anharmonic oscillator. 



where the new coefficients c n (r) are obtained from the old Bender- Wu coefficients c° through 

c n (r) =Ecr( (1 ^] /2 )(2rOr- . 



(10) 



Recall here that the trial frequency f2 was introduced solely as an auxiliary parameter for computational convenience. 
It is not in the original anharmonic potential ([T]), therefore it should not appear in the final result. The dependence 
of Eq. @ can be eliminated by requiring the principle of minimal sensitivity at each order N 
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The VPT ground-state energy is plotted in Fig. |H The axes are the same as in Fig. |3] The black curves are the 
VPT results up to the first four odd orders in the new expansion. They arc almost completely overlapped with each 
other, and one cannot distinguish the difference between different orders by eyesQ However despite of the apparent 
convergence, we still have to ask whether the VPT results converge to the correct value. By performing a strong- 
coupling expansion, It has been shown by Janke and Kleinert [13] that the VPT results for the the first strong-co uplin g 
expansion coefficient ao agreed to all 23 digits with the most accurate value for ao available in the literat ure 
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The convergence radius of the VPT strong-coupling expansion was rigorously proven to be infinity in Ref. |l!6| . 

We have been convinced by VPT in improving the convergence of the weak-coupling expansion from the above 
simple example. The natural question here to ask is whether such trick could be applied to field theory, especially 
to QCD for phenomenological interest. The answer is yes, but we have to do it very carefully, cause simply adding 
a mass term to the QCD Lagrangian would violate gauge symmetry even at the Lagrangian level. We are going to 
address this subtlety in the next section. 



IV. HARD-THERMAL-LOOP PERTURBATION THEORY 



We have seen the oscillating behavior of weak-coupling expansion for thermal gauge theories in Sec. HI In order to 
improve the convergence of perturbation calculations and with the inspiration from VPT, hard-thermal-loop pertur- 
bation theory was introduced by Andersen, Braaten and Strickland as a reorganization of thermal g aug e theory in 
1999 and the one-loop or leading order (LO) thermodynamic calculations were carried out in Refs. [30l432l |. Later 
on HTLpt got extended to two loops or next-to-leading order (NLO) in which a further mass expansion was intro- 
duced to make the calculation tractable in practice [H|,|34j]. The three- loop or next-to-next-to- leading order (NNLO) 
HTLpt calculations have been accomplished recently for both Abclian [35| and non-Abelian [36rl39| theories, and the 
NNLO results turned out to be completely analytic. In the rest of this section, we are going to discuss the setup of 
HTLpt, and results of thermodynamic calculations through NNLO. Please check the above mentioned references for 
calculational details. 



4 The figure is generated by Mathematica for illustration. Due to numerical subtleties, the results up to the first two even orders are not 
easily obtained in Mathematica, therefore the even orders are skipped in Fig. [4] 
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A. Formalism 



In the following of this subsection we will use QCD as an example to show the setup of HTLpt, however the result 
is general for both Abelian and non-Abelian gauge theories as we will comment in the next subsections. 
The Lagrangian density for QCD in Minkowski space reads 

£qcd = -^Tr [G M „GH + ijrfDrf + £ gf + £ gh + A£ QCD , (12) 

where the gluon field strength is G^ v = d^A v — d v — ig[A^, A v \ the term with the quark fields i/j contains an 
implicit sum over the Nf quark flavors, and the covariant derivative is = — igA 1 . The ghost term £ g h 
depends on the gauge-fixing term £ g f. The perturbative expansion in powers of g generates ultraviolet divergences. 
The renormalizability of perturbative QCD guarantees that all divergences in physical quantities can be removed by 
renormalization of the coupling constant a s = g 2 /(4ir) and the necessary counterterms are represented by A£qcd in 
the Lagrangian (112p . There is no need for wavefunction renormalization, because physical quantities arc independent 
of the normalization of the field. There is also no need for renormalization of the gauge parameter, because physical 
quantities are independent of the gauge parameter. 

HTLpt is a reorganization of the perturbation series for thermal gauge theories with the Lagrangian density written 

as 



C= (£qcd + £htl) _ +A£htl. (13) 

g~>VSg 



The HTL improvement term is 




where y M = (l,y) is a light-like four-vector, and (. . ) y represents the average over the directions of y. The term ([14]) 
has the form of the effective Lagrangian that would be induced by a rotationally invariant ensemble of charged sources 
with infinitely high momentum and mod ifies the propagators and vertices self-consistently so that the reorganization is 
manifestly gauge invariant (55l . lll7l4l21 |. The parameter mjj can be identified with the Debye screening mass, and m q 



with the thermal quark mass to account for the screening effects. HTLpt is defined by treating 5 as a formal expansion 
parameter. By coupling the HTL improvement term (|14p to the QCD Lagrangian (fl2|) . HTLpt systematically shifts 
the perturbative expansion from being around an ideal gas of masslcss particles which is the physical picture of the 
weak-coupling expansion, to being around a gas of massive quasiparticlcs which are the more appropriate physical 
degrees of freedom at high temperature. 

Physical observables are calculated in HTLpt by expanding them in powers of (5, truncating at some specified order, 
and then setting 5=1. This defines a reorganization of the perturbation series in which the effects of m D and m 2 
terms in (j!4[) are included to all orders but then systematically subtracted out at higher orders in perturbation theory 
by the 5m 2 D and 5m 2 terms in (|T4l . which is in the spirit of VPT. If we set 5 = 1, the HTLpt Lagrangian (|13p reduces 
to the QCD Lagrangian (fT2j) . If the expansion in 6 could be calculated to all orders the final result would not depend 
on ttid and m q when we set 5=1. However, any truncation of the expansion in 5 produces results that depend on 
Tji y) and TYiq . Some prescription is required to determine mn and m q as a function of T and a s . We will discuss 
several prescriptions in the next subsections. 

The HTL perturbation expansion generates ultraviolet divergences. In QCD perturbation theory, renormalizability 
constrains the ultraviolet divergences to have a form that can be cancelled by the counterterm Lagrangian A£qcd- 
Although the renormalizability of the HTL perturbation expansion has not been proven, the renormalization can 
be archived through NNLO by using only a vacuum counterterm, a Debye mass counterterm and a fermion mass 
counterterm, as well as a coupling constant counterterm. The necessary counterterms for the renormalization through 
NNLO as just discussed read 

A ^o= (T^r + 0(6a s ))(l-5) 2 m 4 D , (15) 



n 128tt 2 £ 

Ato|, = (- UCA 12 J SF a s 5 + 0{5 2 a 2 s )} (1 - S)m% , (16) 

Am^ = f^-±-J±a s 8 + 0{5 2 al)^ (1 - 5)m 2 q , (17) 

6Aa s = - Uc f ~ 4SF a 2 J 2 + 0(S 3 a 3 s ) , (18) 

127T6 



where the coupling constant counterterm is the standard one-loop running of QCD pail2^. 
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B. Thermodynamic potentials through NNLO 



The calculation of the free energy in HTLpt involves the evaluation of vacuum diagrams. In Refs. }30h39I |. the free 
energy was reduced to scalar sum-integrals. The one-loop free energy were evaluated exactly by replacing the sums 
by contour integrals, extracting the poles in e, and then reducing the momentum integrals to integrals that were at 
most two-dimensional and could therefore be easily evaluated numerically. Evaluating two-loop free energy exactly 
would involve the evaluation of five-dimensional numerical integrals which turned out to be intractable, and therefore 
attacking the third loop in this way is hopeless. The fact that mn/q ~ gT suggests that mn /q/ T can be treated as 
expansion parameters of order g in terms of which the sum- integrals can be further expanded [53j . It was shown that 
the first few terms in the m D i q /T expansion of the sum-integrals gave a surprisingly accurate approximation to the 
exact result [30H321 . l53l] . This trick is adopted in the HTLpt results that we are going to present next. 



The Feynman diagrams through NNLO in HTLpt are gathered in Appendix |A"1 Fig. [TU] contains the notation key, 
as well as the LO and NLO HTLpt vacuum diagrams, while Fig. HT1 shows the NNLO HTLpt vacuum diagrams. Note 
that all the propagators and vertices in Figs. [TU] and [TT] arc HTL resummed, which is in contrast to the weak-coupling 
expansion case. The strategy of evaluating the diagrams is to first reduce the diagrams to scalar integrals and then 
expand the resulting scalar integrals in powers of mjj/T and m q /T as discussed above. We will carry out the mrj/T 
and rriq/T expansions to high enough order to include all terms through order g 5 if tod/T and m q /T are taken to be 
of order g. The two-loop approximation will be in the weak-coupling limit accurate to order g 3 and the three-loop 
approximation accurate to order g 5 . We will now present the resulting thermodynamic potentials through NNLO and 



the details of the calculations are presented in Refs. [33H351 1371 . |3_ 



The LO thermodynamic potential reads 



n 



LO 



^idcal 



1 d F 
~Ad~ A 



15 

1 + T ~~ T ^1 



dA 



— 3Q-^m q 



30m 3 



D 



^ log*- 7 - 
4 I 8 2 2 



IE + — 



-4 cn ^ / 2 
m D — 60 — (7T 

dA 



Q)mi , (19) 



where J-ideai is the free energy of a gas of dA masslcss spin-one bosons and tod, rh q and (i are dimcnsionless variables: 



ideal 



m D 



'fq 



m D 



2ttT ' 

TOo 



2ttT ' 
2ttT ■ 



(20) 

(21) 
(22) 
(23) 



The NLO thermodynamic potential reads 



NLO 



•7"id 



= 1 



eal 



1^1 _ 15^3 

Ad A ° 



3tt 

TT 



15 

T 



45to 



45 
T 

D - 



l0£ 



165 / A 
— l0g 2 



IE ~ 

36 

~ IT 



7T" 

"3 



^ + 60^(7r 2 -6)m 
dA 



\ogm D - 2.001 m 2 D + 



-45 (log | + 2.19581J m\ + 180to 



Dm 2 a 



495 



loe 



+ 22 + 7B ) m D 



~ + l5m D + 5 flog | - 2.33452^) m 2 D - 30 flog | - i 



IE 



2 log 2 ) m 3 D 



(24) 
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The NNLO thermodynamic potential reads 



n 



NNLO 



7dt 



15 



1 + --, rtfl 



deal 



4d A 

7T 

fc A a s \' 



D 



3tt 



15 45 „ 
- T + Y m D 



135 



495 



A 5 
l0g 2 + 22 



IE m D 



25 



— m D + 15 I log- 



7b + 2 log 2 I m£, - 90rh z q rh D 



■■in 

1485 
4 

235 A A 
16~ I g 2 



log 



45 1 

4 TO_D 

A 79 

2 ~ 44 



165 / jU 
log g 



7b + log 2 



72 , „ 84 
^logm^-- 



IT 



144 24 
— \ gm D -- lE 



319 
940 4 



315 / u 8 , „ 9 



+ 

111 

235 

m a 

90^- 

TO_D 



3tt 
log 2 



6 74C'(-1) 19 C' (-3) 
Ti lE ~ IT c(-i) + IT C(-3) 



SFOt s \ 15 1 
nil rho 
74 C'(-1) 1 C'(-3) 
47 C(-l) 47 C(-3) 



5 1 25/, fi 1 3 66, „ 4C'(-1) 2 C'(-3) 
4^ + ^r g 2 + 2^ + ^-2^ l0g2+ 5CRy-5CF3T 



7_e + 2 log 2 m_o 



30-r^- 



S2F 



15 45 
_(35 - 321og2)- y m D 



(25) 



Note that the NNLO thermodynamic potential is completely analytic, and the coupling constant counterterm listed 
in Eq. (|18|) coincides with the known one-loop running of the QCD coupling constant 



d/i 



(Uca ~ 4a F )g 4 
24tt 2 



(26) 



In the next subsections we will present thermodynamic quantities as a function of g evaluated at the renormalization 
scale 2nT. 



C. QED thermodynamics 



The thermodynamic potentials for QED are obtained by setting the group factors in Eqs. p^|) . (|2"4")l , and (|2"5)) to 



N f , c A = 



sf = Nf , s 2F = iV/ . (QED) 



(27) 



1. Mass prescriptions 



The mass parameters mo and m/ in HTLpt are in principle completely arbitrary, however physically they are 
proportional to eT at LO in the weak-coupling expansion which has to be satisfied for any mass prescriptions. To 
complete a calculation, it is necessary to specify mo and m/ as functions of e and T. Two possible mass prescriptions 
are considered in Ref. S3 : 



• The variational thermal masses obtained from solving the gap equations 

d 



dmry 
d 

dirif 



J7(T, a,mr),mf,ij,,5 = 1) = 
f2(T, a, mu, m/, [i, 6 = 1) = 



• The e 5 perturbative Dcbyc mass [13, Il24l | and the e 3 perturbative fermion mass [l2 



ml = -N f e 2 T 2 



-N f e 2 T 2 



^2 (47 + 7 + 41og^ + 81og2 



eV3 

47T 3 



2.854 
1 - —a — e 

47T 



(28) 
(29) 

(30) 
(31) 
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FIG. 5. A comparison of the renormalization scale variations between NLO and NNLO HTLpt predictions for the scaled 
pressure of QED with Nf = 1 and the variational thermal masses (left panel) and the perturbative thermal masses (right 
panel). The bands correspond to varying the renormalization scale /z by a factor of 2 around /i = 2nT. 



2. Pressure 



The resulting predictions for the pressure normalized to that of an ideal gas as a function of e evaluated at the 
renormalization scale fj, = 2irT are shown in Fig. [5] Note that when the pressure is evaluated at a scale different than 
fi = 2nT, we use one-loop running of e to determine the value of the coupling at /i = 2irT. As can be seen from 
these figures both the variational and perturbative mass prescriptions seem to be consistent when going from NLO to 
NNLO. At the central value \i = 2irT, both prescriptions arc the same to an accuracy of 0.6% at e = 2.4. As a further 
check, we show a comparison of the NNLO HTLp t res ults with a three-loop calculation obtained previously using 
a truncated three-loop ^-derivable approximation [l26| in Fig. [51 As can be seen from this figure, the agreement 
between the NNLO ^-derivable and HTLpt approaches is at the subperccntage level even at large coupling. The 
improved convergence for QED gives us confidence to apply the same HTLpt reorganization scheme to QCD. 




0.5 1 1.5 2 

e(2jtT) 



FIG. 6. A co mparison of the predictions for the scaled pressure of QED with Nf = 1 between three-loop "^-derivable approxi- 
mation [lH and NNLO HTLpt at n = 2tvT. 
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D. QCD thermodynamics: Pure-glue and Nf = 3 

The thermodynamic potentials for QCD are obtained by setting the group factors in Eqs. and (f2"5j) to 

d A = N 2 C - 1 , d F = N c Nf , c A = N c , s F = ^ , s 2F = ^f^N f . (QCD) (32) 

The pure-glue QCD results are therefore obtained in the Nf = limit 

d,A = N% — 1 , dp = , = 2V C , s F = , S2F = . (pure-glue) (33) 

1. Mass prescriptions 

As pointed out earlier, the HTL mass parameters are completely arbitrary and we need a prescription for them 
in order to complete a calculation. The variational mass prescription unfortunately gives rise to a complex Debye 
mass and m q = at NNLO. One strategy is therefore to throw away the imaginary part of the thermodynamic 
potential to obtain thermodynamic functions that are real valued. The problem of a complex mass parameter was 
encountered at NNLO in SPT massless </> 4 theory [52j], as well as in HTLpt QED [35|, so this seems to be a general 
issue with SPT/HTLpt at NNLO, which is not specific to QCD. The full variational pressure for QED contains both 
the real part as shown in the left panel of Fig. [5] and a imaginary part whose absolute value is about 1% of the real 
part at large coupling. Since the imaginary part is tiny, therefore the real part provides a good approximation to 
the pressure as confirmed by the comparison in Fig. [5] However for QCD, the absolute value of the imaginary part 
of the variational Debye mass becomes comparable to the real part at intermediate coupling, the variational mass 
prescription is therefore not preferab le. Since the weak-co upling per turbative Debye mass receives contributions from 
the nonperturbative magnetic scale |l27l . ll28l | beyond LO |l29l . Il3dl |. for the NNLO HTLpt QCD results that we are 
going to present, the Braaten and Nieto's (BN) mass parameter of three-dimensional electrostatic QCD (EQCD) will 
be used as a substitute for the perturbative Debye mass, effectively discarding the nonperturbative contributions. In 
Ref. [l|l[2(| it was calculated to NLO, giving 



— — T <c A + s F + — — - + — je H log — 1 



3?r V 4 2 2 2 / 1" 

s%a s fl 4 2 2 , A\ 3 s 2F a 



■ ^ ~ 3 log 2 + \ ( 7B + log f 

The perturbative quark mass does not suffer from the nonperturbative magnetic scale, and it turns out that the 
final NNLO results are very insensitive to whether one chooses a perturbative mass prescription for m q , or uses the 
variational mass m q = 0. We will therefore use m q = for simplicity in the following whenever full QCD is concerned. 



2. Pressure 



In the left panel of Fig. [71 we show the LO, NLO, and NNLO HTLpt predictions for the pressure of pure-glue 
QCD normalized to that of an ideal gas as a function of T/T c using the BN mass prescription (f3"4"]l and three -loop 



running of a s [2J]. The points are lattice data for pure-glue with N c = 3 from the Bielefeld collaboration |131 
The bands are obtained by varying the renormalization scale by a factor of 2 around the central value /i = 2irT. 
From the plot we see that the convergence of the successive approximations to the pressure is improved over naive 
perturbation theory. For example, using the naive perturbative approach and comparing the full variation in both 
successive truncations and renormalization scale variation, one finds that at T = 3 T c there is variation in the pressure 
of 0.69 < 'P/'Pidoai < 1-32 Hl|, whereas using HTLpt there is only a variation of 0.74 < V /V ideal < 0.95. Additionally, 
at NNLO we see that the /j, = 2nT result for the pressure in Fig. [7J coincides with the lattice data down to T ~ 3T C . 

In the right panel of Fig. [Jj we show the LO, NLO, and NNLO HTLpt predictions for the pressure of QCD 
with Nf = 3 normalized to that of an ideal gas as a function of T using the BN mass (fM]) as well as m q = 0. 
For the strong coupling constant a s , we use three- loop running [24| with A^jg = 344 MeV which for Nf = 3 gives 
a s (5GeV) = 0.2034 [23. The band is again obtained by varying the renormalization scale by a factor of 2 around 
the central value /z = 2irT. The Nf = 2 + 1 lattice data from the Wuppertal-Budapest collaboration use the stout 
action. Since their results show essentially no dependence on the lattice spacing (it is smaller than the statistical 
errors) , they provide a continuum estimate by averaging the trace anomaly measured using their two smallest lattice 
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FIG. 7. Comparison of LO , NL O, and NNLO predictions for the scaled pressure for pure-glue QCD with lattice data from 
the Bielefeld collaboration |13ll | (left panel) and Nf = 3 QCD with Nf = 2 + 1 lattice data from the hotQCD [l32| and 
Wuppertal-Budapest [l33| collaborations (right panel). Shaded bands show the result of varying the renormalization scale \x 
by a factor of 2 around /i — 2nT for the NNLO result. 



spac ings corresponding to N T = 8 and N T = 10 [133l |. which were essentially on top of the N T = 6 measurement 
jl34l ]Pl Using standard lattice techniques, the continuum-estimated pressure is computed from an integral of the trace 
anomaly. The Nf = 2 + 1 lattice data from the hotQCD collaboration are their N T = 8 results using both the asqtad 
and p4 actions [132| . The hotQCD results have not been continuum extrapolated and the error bars correspond to 
only statistical errors and do not factor in the systematic error associated with the calculation which, for the pressure, 
is estimated by the hotQCD collaboration to be between 5 - 10%. As can be seen from the right panel of Fig. 
the successive HTLpt approximations represent an improvement over the successive approximations coming from a 
weak-coupling expansion; however, as in the pure-glue case in the left panel of Fig. [7J the NNLO result represents a 
significant correction to the LO and NLO results. That being said, the NNLO HTLpt result agrees quite well with 
the available lattice data down to temperatures on the order of 2T C ~ 340 MeV for QCD with Nf = 30 Below these 
temperatures the successive approximations give large corrections with the correction from NLO to NNLO reaching 
100% near T c . 



3. T 4 scaled trace anomaly 

In the left panel of Fig. [5] we show the NNLO HTLpt prediction for the trace anomaly of pure-gluc QCD normalized 
to T 4 as a fu nctio n of T. The points are lattice data for pure-glue with N c = 3 from the Wuppertal-Budapest 
collaboration |l35| . For temperatures below approximately 2 T c ~ 500 MeV, there is a large discrepancy between the 
HTLpt prediction and lattice data. The discrepancy decreases as increasing temperature, and for temperatures above 
approximately 4T C ~ 1000 MeV, the NNLO HTLpt result is in good agreement with the lattice result. 

In the right panel of Fig. [5] we show the NNLO HTLpt prediction for the trace anomaly of QCD with Nf = 3 
normalized to T 4 as a function of T. The data from both the Wuppertal-Budapest collaboration and the hotQCD 
collaboration are taken from the same data sets displayed in the right panel of Fig. [7J and described previously. In 
the case of the hotQCD results we note that the results for the trace anomaly using the p4 action show large lattice 
size affects at all temperatures shown and the asqtad results for the trace anomaly show large lattice size effects for 
T > 200 MeV. We see very good agreement between the HTLpt prediction and the available lattice data down to 
temperatures on the order of T ~ 2 T c . 



5 It is noted that the Wuppertal-Budapest group has published a few data points for the trace anomaly with N T = 12 and within statistical 
error bars these are consistent with the published continuum estimated results. 

6 The Wuppertal-Budapest and hotQCD data were obtained using a physical strange quark mass; however, HTLpt calculations use 
masslcss quarks. The difference between massive and masslcss quarks is expected to be significant only for T < 32 MeV corresponding 
to the temperature where the lowest fermionic Matsubara mode equals the strange quark mass. 
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FIG. 8. Comparison of NNLO pr edict ions for the T 4 scaled trace anomaly for pure-glue QCD with lattice data fro m the 
Wuppertal-Budapest coll abora tion |l35l ] (left panel) and Nf = 3 QCD with Nf — 2 + 1 lattice data from the hotQCD [l32t | 
and Wuppertal-Budapest [l33f | collaborations (right panel). Shaded band shows the result of varying the renormalization scale 
fj, by a factor of 2 around fi = 2nT for the NNLO result. 



4- T scaled trace anomaly 

In the left panel of Fig. |H1 we show comparison of the trace anomal y sca led by T 2 T 2 for pure-glue QCD bet ween 
NNLO HTLpt results and lattice data from the Wuppertal-Budapest pi], Bielefeld |l3l| , and WHOT-QCD [l36j 
collaborations. The solid black line is the NNLO HTLpt result obtained using a one-loop running coupling and the 
dashed black line is the HTLpt result obtained using a three-loop running coup ling. In the case of the three-loop 
running the lattice determination of T c /Aj^g = 1.26 is used to fix the scale [135j . For comparison between the one- 
and three-loop results we require that both give the same value for the strong coupling when the renormalization 
scale /j, = 5 GeV. Numerically, one finds a s (5 GeV) = 0.140553. The difference between one- and three-loop running 
will be used to gauge the theoretical uncertainty of the NNLO HTLpt results. For both the one- and three-loop 
running /i = 2irT is taken. As can be seen from the plot, at low temperatures NNLO HTLpt underpredicts the trace 
anomaly for pure-glue QCD and only starts to agree at temperatures on the order of 8T C . At high temperatures 
one sees excellent agreement between the NNLO HTLpt results and the lattice results. We note in this context 
that if one allows for a fit to the unknown perturbative g 6 contribution to the pressure, then the EQCD approach 
can equally- well describe the pure Yang- Mills lattice data down to temperatures on the order of 8 T c [l37l . Il38l | . 
The most remarkable feature of the lattice data for the T 2 scaled trace anomaly is that it is essentially constant in 
the temperature range T ~ 1.1 — 4T C . It has been suggested that this behavior is due to the influ ence of p ower 
corrections of order T 2 which are nonperturbative in nature and might be related to confinement [139l4l44| . At 
temperatures above approxim ately 4 T c the latest Wuppertal-Budapest results show an upward trend in accordance 
with perturbative predictions [135| . The WHOT-QCD results also exhibit an upward trend, however, it starts at 
much lower temperatures. This discrepancy coul d be due to their fixed scale approach not having sufficiently large 



N T at high temperatures as noted in their paper 136 1. 



In the right panel of Fig. [9j we show the NNLO HTLpt trace anomaly scale d by T 2 for Q CD with N / = 3 and 
compare to the Nj = 2 + 1 lattice results available from the Wuppertal-Budapest [133| . hotQCD [I32ill45j . and RBC- 
Biclcfeld jl46lll47| collaborations. As before, we show HTLpt results using both one- and three-loop running couplings 
with the requirement that both couplings give Qs(5 GeV) = 0.2034 in accordance with recent high precision lattice 
measurements of the running coupling constant [25|. The lattice data from the Wuppertal-Budapest collaboration 
are taken from the same data sets displayed in Figs. [7] and [5] and described previously. The lattice data from the 
hotQCD colla boration are their N T = 8 results using the asqtad, p4, and HISQ actions which have not been continuum 
extrapolated jl32l . 1145] . The latti ce data from the RBC-Bielefeld collaboration is N T = 6 and have also not been 
continuum extrapolated jl46l . Il47j . As can be seen from the plot for T > 400 MeV, one finds reasonable agreement 
between the NNLO HTLpt predictions and available lattice results. Below this temperature the Wuppertal-Budapest 
and hotQCD results do not seem to agree. Therefore, it is difficult to draw conclusions about the efficacy of the HTLpt 
approach; however, naively one expects that for temperatures less than twice the critical temperature T c (^170MeV 
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FIG. 9. Comparison of N NLO pre dictio ns for the T 2 scale d tr ace anomaly for pure-glue QCD with lattice data from the 
Wuppertal-Budapest [HI], Bielefeld [HH, and WHOT-Q CD [l3l c ollaborations (left panel) and N f = 3 QCD with N f = 2 + 1 
lattice data from the Wuppertal-Budapest [l33T ]. hotQCD [l32l . ll45l ]. and RBC-Bielefeld [Tia . [l47l ] collaborations (right panel). 
The HTLpt results are taken at /i = 2nT. 



for QCD with Nf = 2 + 1) that nonperturbative corrections should become increasingly important. 



5. Discussions 



The HTLpt results indicate that the lattice data at temperatures T ~ 2 T c are consistent with the quasiparticle pic- 
ture. This is a nontrivial result since, in this temperature regime, the QCD coupling constant is neither infinitesimally 
weak nor infinitely strong with g ~ 2, or equivalently a s ~ 0.4. Therefore, we have a crucial test of the quasiparticle 
picture in the intermediate coupling regime. 

The failure of HTLpt to match lattice data at lower temperatures is to be expected since one is expanding around 



the trivial vacuum = and therefore neglec ts the app roximate center symmetry Z(N C ), which becomes essential 
as one approaches the deconfinement transition |l48l - ll54T |. In addition, it is also in line with expectations since below 
T ~ 2 — 3 T c a s imple "electric" quasiparticle approximation breaks down due to nonpertu rbative chromomagnetic 



onpe 

effects |l27lll28j . Besid es, t here have been also hints t hat gauge- fixing ambiguities [l55HT"6l| . topological objects such 
as quantum instantons [162| and magnetic monopoles |163l Il64j might play important roles on the thermodynamics 
at intermediate temperature. 

We find that when including quarks the agreement with lattice data is greatly improved as compared to the NNLO 
results of pure-glue QCD. Fermions are perturbative in the sense that they decouple in the dimensional-reduction 
step of effective field theory, so we expect that including contributions from quarks gives at least as good agreement 
with the lattice calculations as the pure-glue case. However, the exact reason for the better agreement between the 
HTLpt predictions and lattice calculations when including quarks is not clear to us. 

Just as for NNLO QED and massless scalar c/> 4 theory, we encountered again the complex variational Debye mass 
when solving the gap equations. Whether the complexity of the variational Debye mass parameter is due to the 
additional expansion in mp/T and m q /T is impossible to decide at this stage. The correction to the pressure going 
from NLO to NNLO is also rather large. It is unfortunate that the nonperturbative magnetic scale prevents going to 
N 3 LO without supplementing the calculation with input from three-dimensional lattice calculations, as it would be 
interesting to see whether the complexity of the Debye mass parameter and the slow convergence persists. 



V. CONCLUSIONS AND OUTLOOK 



In this paper, wc have briefly reviewed the progress of hard-thermal-loop perturbation theory that has been made 
over the past 12 years concentrating on systematic thermodynamic calculations. We began by using the quantum 



15 



mechanical anharmonic oscillator to show the breakdown of weak-coupling expansion which demonstrated, in sharp 
contrast to our intuition, that small coupling expansion may not be the same as perturbative expansion, or small 
coupling 7^ perturbative. Then we showed that as a rcsummation scheme, VPT was able to convert divergent weak- 
coupling expansions into fast converging approximations. With the inspiration from VPT, we introduced HTLpt as a 
reorganization for thermal gauge theory. The thermodynamic application of HTLpt leads to laudable results for both 
Abelian and non-Abelian theories. 

The success of HTLpt is not totally unexpected since it is essentially just a reorganization of perturbation theory 
which shifts the expansion from around an ideal gas of masslcss particles to that of massive quasiparticlcs which 
are the more appropriate degrees of freedom at high temperature. The fact that the mass parameters are not 
arbitrary but functions of g and T determined variationally or perturbatively also indicates that HTLpt does not 
modify the original gauge theory but just reorganizes its perturbation series. Gauge invariance which is guaranteed 
by construction in HTLpt is useful both as a consistency check in calculations and as a way to simplify calculations. 
Although the renormalizability of HTLpt is not yet proven, the fact that it is renormalizable at NNLO using only 
known countcrterms shows promising light along the way. 

So far, thermodynamics for quantum fields has been studied intensively in the community, both perturbatively 
through higher orders or numerically on the lattice, however real-time dynamics is still in its very early stage of 
development. Transport coefficients arc of great interest since they are theoretically clean and well defined non- 
equilibrium dynamical quantities. Along the line of pert urbative approach to transport coefficients, the only known 
ones to NLO are s hear visco sity in scalar <fi 4 theory [165J, heavy quar k dif fusion in QCD and M — 4 supcrsymmetric 
Yang-Mills theory [l66l Il67j , and transverse diffusion rate q in QCD jl68l | , and all of them exhibit poor convergence 
as bad as the case of thermodynamic quantities, such as the pressure. Since dynamical quantities are still not well 
described by lattice gauge theory, new resummation techniques are urgently needed in order to achieve a better 
understanding of transport coefficients. 

Although the papers written to date have focussed on using HTLpt to compute thermodynamic observables, the goal 
of this project is to create a framework which can be applied to both equilibrium and non-equilibrium systems. HTLpt 
is formulated in Minkowski space, so its application to real-time dynamics is straightforward. With the confidence 
from thermodynamic calculations, HTLpt is ready to enter the domain of real-time dynamics at temperatures that 
are relevant for LHC and this might be of great help in deepening our knowledge in the properties of the quark-gluon 
plasma. Last but not least, it would be also interesting to explor e the applicability of the HTLpt /SP T method to 
other finite temperature or density systems, such as compact stars jl69j | and ultracold atoms [l7Ct Il7l| . 
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FIG. 10. QCD diagrams contributing through NbO in HTLpt. The spiral lines are gluon propagators, the dotted lines are ghost 
propagators, and the solid lines are quark propagators. A circle with a II indicates a one-loop gluon self-energy insertion, a E 
indicates a one-loop quark self-energy insertion, and a F indicates a one-loop vertex insertion. A square with a g is shorthand 
for the pure-glue diagrams contributing to the one-loop gluon self-energy. All gluon and quark propagators and vertices shown 
are HTL-resummed propagators and vertices. The logic behind the diagram notation is as follows: diagrams consisting only of 
gauge propagators have g superscripts. Diagrams containing fermion propagators have / superscripts. The subscript indices 
are identical to those used in [37| (pure-glue QCD) and [35| (QED). We do not display the symmetry factors in the diagrams. 



